Nuclear spin-lattice relaxation rate in the D+iD superconducting state: implications 

for CoO superconductor 
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We calculated the nuclear spin-lattice relaxation rate 1/Ti for the D+iD superconducting state 
with impurities. We found that small amount of unitary impurities quickly produces the residual 
density of states inside the gap. As a result, the T-linear behavior in 1/Ti is observed at low 
temperatures. Our results show that the D+iD pairing symmetry of the superconducting state of 
Nao.35Co02 ■ 1/H2 O is compatible with recent 59 Co 1/Ti experiments of several groups. 

PACS numbers: 74.20, 74.20-z,74.50 



I. INTRODUCTION 

The recent discovery of Nao.35Co02 • 2/H2 O 
superconductor 1 has renewed interest in the study of ex- 
otic superconductivity (SC). Motivated by the triangu- 
lar lattice structure of spin-i Co atoms of the reference 
C0O2 layer, several authors 2 have quickly proposed that 
it might be a realization of the long sought spin-liquid 
based superconductivity, originally proposed for high-T c 
superconductors by Anderson et al. 3 . Indeed several the- 
oretical calculations 4 indicate that the frustration of the 
antifcrromagnctism (AFM) on the triangular lattice of 
Co favors the D+iD superconducting symmetry. 

For a new superconducting material, the first step to- 
ward laying down the theory of the superconducting state 
is the experimental identification of the SC gap symme- 
try. Nuclear spin-lattice relation rate 1/Ti and Knight 
shift measurements have been the successful tools for this 
purpose and several measurements of 1/Ti and Knight 
shift of this compound have already been reported. At 
present, however, these reports appear inconsistent with 
each other. Regarding the coherence peak of Co 59 1/Ti, 
for example, Kobayashi et al. 5 and Waki et al. 6 reported 
the existence of a coherence peak; the maximum value 
of 1/TiT is about 1.1 to 1.5 times of the normal state 
value while an ideal s-wave case has almost twice of the 
normal state value 7 . However, Fujimoto et al. 8 recently 
reported no coherence peak. Furthermore, they reported 
more precise 59 Co 1 /Ti below T c and found a substantial 
region of linear temperature dependence in 1 /Ti at low 
temperatures and T 3 behavior near and below T c . To fit 
their data a residual density of states (DOS) of 0.65N(0) 
was needed assuming an impure D-wave SC state. From 
this calculation, these authors concluded that the non S- 
wave gap with the lines of nodes is most consistent with 
their 1/Ti data and the compatibility of the D+iD sym- 
metry gap having a full gap remains as a question. 

In this note, we calculated 1/Ti for D+iD gap with 
impurities using self-consistent T-matrix approximation 
and quantitatively compared our results with experi- 
ments reported up to date. Since the impurity effects 
on various symmetry gaps has been studied by many 



authors 9,10 and the generic behaviors are well known, we 
summarize them. Within T-matrix approximation, the 
impurity scattering renormalizes the self-energies both of 
normal and anomalous components. In Nambu notation, 
they correspond to £ , £2, and £3, respectively; £0 an d 
£3 are the two separate components of the normal self- 
energy according to particle-hole transformation and £2 
is the anomalous self-energy component. Among them, 
the renormalization of £3 usually vanishes for particle- 
hole symmetric band, which is a good approximation for 
most cases. For an isotropic s-wave gap, both £0 and 
£2 are renormalized and their effects largely cancel each 
other in the renormalized gap equation; as a result the 
impurity effect is null for low temperature SC properties 
such as the changes of the order parameter (OP) magni- 
tude (Ao) and T c , and the residual DOS etc. For a sign 
changing OP such as the D-wave gap, the renormaliza- 
tion of £2 vanishes, too, due to symmetry, and the impu- 
rity effects appear only in the renormalization of £0 the 
normal self-energy. The case of D+iD gap is interesting 
in that it has a full gap all over the FS and should dis- 
play similar thermodynamic behaviors as in an isotropic 
s-wavc gap. However, regarding the impurity effects, it 
is expected to show a mixed behavior of a D-wave gap 
and a s-wave gap because the sign changing OP leads to 
a vanishing £2 renormalization despite a full gap. 

We find that for unitary impurities the residual DOS is 
created in the D+iD gap as easily as in the D-wave gap. 
This happens because the resonant scattering process of 
unitary impurity is more efficient at lower DOS. As a re- 
sult, D+iD gap displays a similar T-linear behavior in 
1 /Ti at low temperatures with the same amount of uni- 
tary impurities as in the simple D-wave case. Therefore, 
the observation of the T-linear behavior in 1/Ti at low 
temperatures 8 itself cannot rule out the D+iD symmetry 
gap on C0O2 superconductor. Moreover, in contrast to 
the D-wave gap the D+iD gap with a smaller amount 
of impurities displays a reduced coherence peak reflect- 
ing full opening of the gap around the FS. Putting to- 
gether our theoretical results, we conclude that the D+iD 
symmetry gap with unitary impurities provides the most 
consistent agreement with the currently available nuclear 
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spin-lattice relaxation data. We also studied the effects 
of Born limit impurity on both symmetry gaps. We find 
that the response to Born limit impurities is qualitatively 
different for D+iD and D-wave symmetry gaps, and both 
cases cannot fit the experimental data. 



II. FORMALISM 

The effect of the impurity scattering is included with 
T-matrix approximation 9 . As we explained above, for 
sign-changing OP and assuming particle-hole symmet- 
ric band, we just need to calculate So = TTo, where 
r = rii/irNo, No the normal DOS at the Fermi energy, 
rii the impurity concentration. To is defined in Matsub- 
ara frequency as 



To(ui n ) 



(i) 



where Co = lo + T,q(lo) and (...) means the average over 
the FS. The first term in the bracket of Eq. (4) is N 2 (lo). 
The second term vanishes in our calculations because 
of the symmetry of the OP. To calculate \jT\T using 
Eq.(4), we need the full temperature dependence of the 
gap function A(k). The gap equation Eq. (3) is basically 
the BCS gap equation, therefore it gives the BCS tem- 
perature dependence for A(k) and Ao/T sc = 1.764 and 
2.14 for the D+iD- and D-wave gap solutions, respec- 
tively. We use the phenomenological formula, A(k,T) = 
A(k,T = 0) E(T); g(T) = tanh(/3 y/T sc /T - 1), with pa- 
rameters /? and Ao/T sc . The temperature dependence 
of So(w,T) (= r To(w,T)) is similarly extrapolated: 
T (w,T) = T (lo,0) S(T) + T normal (l - 3(T)), where 
Tnormai = r '/ 1 (c 2 + 1) is the normal state To. Then we 
only need to calculate A(k) and Tq at zero temperature. 



III. RESULTS AND DISCUSSIONS 



and 



nN ^ Col 



A 2 (k)' 



where Co n = to n + So and the scattering strength param- 
eter c is related to the s-wave phase shift 5 as c = cot(<5). 
With this To the following gap equation is solved self- 
consistently. 



A(k) = -N J2v(k,k') 



de 



A(k 



D Q 2 l + e 2 + A 2 (k')' 



For D+iD gap equation, we assume the pairing potential 
V(k, k') to be a constant because thermodynamic prop- 
erties of D+iD SC are determined by the absolute magni- 
tude of the gap \D+iD\ which is the same as the isotropic 
s-wavc gap. The property of the D+iD gap in the above 
gap equation manifests itself only as the absence of £2 
renormalization. For comparison, we also calculate the 
D-wave case with a proper pairing potential 11 , which also 
has no £2 renormalization due to the sign-changing OP. 
The only technical difference of the impurity effects on 
D+iD- and D-wave symmetry gaps is the normal self- 
energy correction due to impurities in the SC states with 
a full gap and lines of nodes in the gap, respectively. 

With the gap function A(k) and Tq(lo) obtained from 
Eq.(l) (Tq(lo) is analytically continued from To(w„) by 
Pade approximant method), we calculate the 1/Ti nu- 
clear spin-lattice relaxation rate 9,12 
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gMg) 
dw 



((Re 



to 



^Cu 2 -A 2 (k) 



((Re - AW )..)'] 
U ^JCo 2 - A 2 (k)' ' J 



(2) In our numerical calculation, the BCS value (3 — 

1.74 and the ratio 2Ao/T sc = 3.5 for D+iD gap and 
2 Ao/T^ = 5 for D-wave gap are used, respectively. Note 
that the final results are insensitive to the precise value 
of p. Fig. 1(a) shows the normalized nuclear spin-lattice 
relaxation rate 1 /Ti for D+iD SC gap with varying impu- 
rity concentration of unitary impurities (c=0). The nor- 
malization is chosen so that the value of 1 /Ti at T c is 100 
in arbitrary units. The experimental data (black circles) 
are also normalized in the same fashion. The inset shows 
the corresponding DOS. The key result we find is that 

(3) 

although D+iD gap develops a full gap as in an isotropic 
s-wave gap, the unitary impurity scattering produces res- 
onance states inside the gap (at lo = in our case of c=0 
and a particle-hole symmetric band) . A finite amount of 
unitary impurities, therefore, quickly produces a resid- 
ual DOS around lo = inside the gap. As a result, the 
nuclear spin-lattice relaxation rate 1 /Ti develops a grow- 
ing T-linear region at low temperatures with increasing 
impurity concentration. The high-temperature 1/Ti be- 
low T c initially follows the pure D+iD curve (black open 
squares) with a small amount of impurities, but with in- 
creasing impurity concentration it follows a nongeneric 
power law dependence 13 . With the impurity concentra- 
tion of r/Ao = 0.32 (Ao is the absolute magnitude of the 
pure D+iD gap at zero temperature), the residual DOS 
at lo = N(0) is about 0.63 No (No is the normal state 
DOS at Fermi level). Compared with the experimental 
data of 59 Co 1/Ti from T. Fujimoto et al. (Ref.[8]) it 
provides a reasonably good fit: the low temperature T- 
linear behavior and its magnitude, the high temperature 
power law (experimental fit was T 2 - 2 ), and the absence 
of the coherence peak. Therefore, the experimental ob- 
servation of the substantial region of T-linear behavior 8 
/ 4 \ in 59 Co 1/Ti of Na .3 5 CoO 2 • J/H 2 O doesn't rule out the 
possibility of D+iD gap in this compound. 



2 



Furthermore, our calculations for D+iD gap provide 
a resolution for the conflicting observation of the coher- 
ence peak by other groups 5,6 . First, without impurities 
(black open squares) the height of the coherence peak is 
about 1.46 times of the value of 1/Ti at T c , which is 
a little smaller than the value for isotropic s-wave case. 
This is because the second term in Eq.(4) vanishes due 
to the symmetry of the D+iD OP. With impurities the 
coherence peak is quickly suppressed and there is a small 
peak visible with a smaller amount of impurities (for ex- 
ample, sec T/Ao = 0.04 (open purple squares); note the 
logarithmic scale). We notice that there is a systematic 
trend of the height of the coherence peak with T c of the 
samples used for 59 Co 1/Ti measurement: no coherence 
peak with T c =3.9 K 8 ; coherence peak with height 1.1 
times of the value of 1/Ti at T c with T c =4.5 K 5 ; co- 
herence peak height 1.4 times of the value of 1/Ti at T c 
with T c =5.0 K 6 . Therefore, we can understand the seem- 
ingly inconsistent experimental reports on 59 Co 1/Ti by 
different groups with D+iD gap and as a result of the dif- 
ferent samples with different scattering rates 14 . We also 
show the D-wave results with unitary impurities (c=0) in 
Fig. 1(b). The overall shape of DOS for D-wave gap with 
and without impurities is qualitatively different from the 
case of D+iD gap. Nevertheless, the accumulation of 
the residual DOS N(0) with impurities arc very similar 
in magnitude, which is the main property determining 
the low temperature 1/Ti behavior. As a result, we can 
also fit the data of Fujimoto et al. with a D-wave gap 
with a similar amount of impurities (the necessary T/Aq 
value should be between 0.32 and 0.16, see Fig. 1(b)); the 
residual DOS N(0) with r/A =0.32 is about 0.73 N . 
However, the D-wave results cannot explain the obser- 
vation of the coherence peak seen by some experimental 
groups 5,6 . For completeness, we show the same calcula- 
tions in the Born limit (c=l) in Fig. 2(a) and Fig. 2(b). 
Both D+iD and D-wave gaps with Born limit impurities 
doesn't produce a close fit to the experimental data un- 
less r/A ->• 0.8. Even with r/A = 0.8, the detailed 
line shape of the calculated 1/Ti doesn't quite fit the 
experimental data for both gaps. Moreover, the reduc- 
tion of T c and the gap magnitude Ao with T/Aq = 0.8 
is about 50 % or more for both gap cases, which is not 
consistent with experiments. 



low temperatures and can fit the experimental data of 
Fujimoto et al. 8 with the unitary impurity concentration 
of T/Ao ~ 0.3. However, considering observations of the 
coherence peak in 59 Co 1/Ti by other groups 5,6 and a 
systematic trend of the height of the coherence peak with 
T c of the samples, we conclude that the D+iD symmetry 
gap is the most consistent with the currently available 
nuclear spin-lattice relaxation data of Nao.3sCo02 ■ j/H 2 
O. 
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In summary, we have calculated the nuclear spin-lattice 
relaxation rate 1/Ti for both D+iD- and D-wave super- 
conducting states with impurities. We found that unitary 
impurities produce the residual density of states inside 
the gap due to the resonant scattering for both cases. 
Interestingly, we found that the rate of accumulation of 
the residual density of states with unitary impurities is 
similar for both SC gaps. As a result, both D+iD and 
D-wave gaps produce the T-linear behavior in 1/Ti at 
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FIG. 1. (A) The normalized 1/Ti of D+iD gap with uni- 
tary impurities (c=0). Experimental data (black circles) are 
also normalized (T. Fujimoto et al. (Ref[8])). Inset is the cor- 
responding DOS. (B) The same plots as in (A) for a D-wave 
gap with unitary impurities (c=0). 



FIG. 2. (A) The normalized 1/Ti of D+iD gap with Born 
limit impurities (c=l). Experimental data (black circles) are 
also normalized (T. Fujimoto et al. (Ref[8])). Inset is the cor- 
responding DOS. (B) The same plots as in (A) for a D-wave 
gap with Born limit impurities (c=l). 
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peratures near T c is not a generic property for any sym- 
metry gaps; it is determined by the temperature slope of 
A(T) near T c , elastic and inelastic scattering rates, and 
the anisotropy of the gap, etc. Often quoted T 3 behav- 
ior of 1/Ti for a pure D-wave gap is a generic property 
at low temperatures reflecting the lines of nodes, but the 
same T 3 behavior observed at high temperatures is a re- 
sult from combined effects mentioned above. Nevertheless, 
empirically once the coherence peak is suppressed the high 
temperature power law of 1/Ti ~ T a can be fit with a 
between 2 and 3 for the dirty D+iD and D cases. 
This conclusion needs further experimental verification 
since the measurements of Ref.[5] and Ref[6] have some 
ambiguity in extracting the Ti relaxation time. 
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